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Abstract 

We study, both classically and quantum-mechanically, the problem 
of a neutral particle with a spin S, mass m and magnetic moment fi, 
moving in one dimension in an inhomogeneous magnetic field given by 

B = B G i+B' ± xy. 

This problem serves for us as a toy model to study the trapping of 
neutral particles. We identify 

which is the ratio between the precessional frequency of the particle 
and its vibration frequency, as the relevant parameter of the problem. 

Classically, we find that when is antiparallel to B, the particle 
is trapped provided that K < 0.5. We also find that viscous friction, 
be it translational or precessional, destabilizes the system. 

Quantum-mechanically, we study the problem of a spin S = h/2 
particle in the same field. Treating K as a small parameter for the 

*Also with the Department of Physics, University of California, San Diego, La Jolla, 
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perturbation from the adiabatic Hamiltonian, we find that the lifetime 
T esc of the particle in its trapped ground-state is 
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cxp 



2K 



where T v ^ = 27ry mBo/fi [B'jj is the classical period of the particle 
when placed in the adiabatic potential V = \i |B|. 



1 Introduction. 



1.1 Magnetic traps for neutral particles. 

Recently there has been rapid progress in techniques for trapping samples 
of neutral atoms at elevated densities and extremely low temperatures. The 
development of magnetic and optical traps for atoms has proceeded in parallel 
in recent years. While optical methods have proved to be an efficient means 
of cooling atoms to temperatures in the microKelvin range, further progress 
is limited by interatomic interactions induced by the scattering of photons. 
The effort to attain higher densities and lower temperatures has therefore 
concentrated on the development of purely magnetic traps|], 0, ||, [|, || ||, |7j . 
Such traps exploit the interaction of the magnetic moment of the atom with 
the inhomogeneous magnetic field to provide spatial confinement. 

Microscopic particles are not the only candidates for magnetic traps. In 
fact, a vivid demonstration of trapping large-scale objects is the hovering 
magnetic top|| H |Hj]. This ingenious magnetic device, which hovers in 
mid-air for about 2 minutes, has been studied recently by several authors 
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1.2 Qualitative description. 

The physical mechanism underlying the operation of magnetic traps is the 
adiabatic principle. The common way to describe their operation is in terms 
of classical mechanics: As the particle is released into the trap, its magnetic 
moment points antiparallel to the direction of the magnetic field. While 
inside the trap, the particle experiences lateral oscillations which are 
slow compared to its precession uj prec . Under this condition the spin of the 
particle may be considered as experiencing a slowly rotating magnetic field. 
Thus, the spin precesses around the local direction of the magnetic field 
B (adiabatic approximation) and, on the average, its magnetic moment /i 
points antiparallel to the local magnetic field lines. Hence, the magnetic 
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energy, which is normally given by —\i ■ B, is now given (for small precession 
angle) by [i |B|. Thus, the overall effective potential seen by the particle is 



In the adiabatic approximation, the spin degree of freedom is rigidly coupled 
to the translational degrees of freedom, and is already incorporated in Eq.(^). 
Thus, under the adiabatic approximation, the particle may be considered as 
having only translational degrees of freedom. When the strength of the mag- 
netic field possesses a minimum, the effective potential becomes attractive 
near that minimum and the whole apparatus acts as a trap. To prevent spin- 
flip (Majorana transitions), most magneto-static traps include a bias field, 
so that the effective potential V e ff possesses a nonvanishing minimum. 

As mentioned above, the adiabatic approximation holds whenever uo prec ^> 
uj V ib- As uJ vr ec is inversely proportional to the spin, this inequality can be 
satisfied provided that the spin of the particle is small enough. If, on the 
other hand, the spin of the particle is too large, it cannot respond fast enough 
to the changes of the direction of the magnetic field. In this limit u) prec ^> u vi b, 
the spin has to be considered as fixed in space and, according to Earnshaw's 
theorem [pj|], becomes unstable against translations. 

1.3 The purpose and structure of this paper. 

The discussion of magnetic traps in the literature is, almost entirely, done 
in terms of classical mechanics. In microscopic systems, however, quantum 
effects become dominant, and in these cases quantum mechanics is suited 
for the description of the trap. An even more interesting issue is the under- 
standing of how the classical and quantum descriptions of a given system are 
related. In this paper we study, both classically and quantum-mechanically, 
the quantitative nature of magnetic traps. In order to keep the underly- 
ing physics transparent, we devise a simplified model for the inhomogeneous 
magnetic field of such traps. We further neglect the effect of interactions 
between the particles in the trap and so we analyze the dynamics of a single 
particle inside the trap. For simplicity, the particle is considered to have only 
a single translational degree of freedom. Its spin degree of freedom, on the 
other hand, is taken completely into account. 

The structure of this paper is as follows: In Sec.0 we start by defining the 
system we study, together with useful parameters that will be used through- 
out this paper. Next, we carry out a classical analysis of the problem in 
Sec.[| Here, we find two stationary solutions for the particle inside the trap. 
One of them corresponds to a state whose spin is parallel to the direction of 
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the magnetic field whereas the other one corresponds to a state whose spin 
is antiparallel to that direction. When considering the dynamical stability 
of these solutions, we find that only the antiparallel stationary solution is 
stable. We also study the same problem but with viscous friction added, and 
arrive at the interesting result that friction destabilizes the system. In Sec.|] 
we reconsider the problem, from a quantum-mechanical point of view. Here, 
we also find states that refer to parallel and antiparallel orientations of the 
spin, one of them being bound while the other one unbounded. In this case, 
however, these two states are coupled due to the inhomogeneity of the field, 
and we move on to calculate the lifetime of the bound state. Finally, in Sec.|5| 
we compare the results of the classical analysis with those of the quantum 
analysis and comment on their implications for practical magnetic traps. 



2 Description of the problem. 

We consider a particle of mass m, magnetic moment /i and intrinsic spin S 
(aligned with n) moving in ID space in an inhomogeneous magnetic field B 
given by 

B =B z+B' ± xy. (2) 

This field possesses a nonzero minimum of amplitude at the origin, which is 
the essential part of the trap. Note also that the direction of the field twists 
( or curls) as one moves along the x axis. The Hamiltonian for this system is 

H = ^-,B (3) 

where P is the momentum of the particle. 

We define uj prec as the precessional frequency of the particle when it is at 
the origin x — 0. Since at that point the magnetic field is B =Bqz we find 
that 

= (A\ 
^prec — g ■ [fx) 

Next, we define u^b as the small-amplitude vibrational frequency of the par- 
ticle when it is placed in the adiabatic potential field given by 

V(x) = \i \B(x) | = fiB + ± 2 x 2 J + O (x 4 ) . 
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For this potential we have 

, d 2 V 



dx 2 

and therefore 



B' ) 2 



x=0 B 



k _H& L f l i 



u vib =\— = \ — . (5) 

m v ttlBq 



We also define the ratio between uj^ and u prec , 



Ks <*» = my, (6) 

This will be our 'measure of adiabaticity'. It is clear that as K becomes 
smaller and smaller, the adiabatic approximation becomes more and more 
accurate. Note also that K is the only possibility to form a non-dimensional 
quantity (up to an arbitrary power) out of the parameters of the system. 
The value of K therefore, completely determines the behavior of the system. 



3 Classical analysis. 

3.1 The stationary solutions. 

We denote by n a unit vector in the direction of the spin (and the magnetic 
moment). Thus, the equation of motion for the center of mass of the particle 
is 

m ^ = /i ^ (n ' B) ' (7) 

and the evolution of its spin is determined by 

S^unx B. (8) 
dt 

The two equilibrium solutions to Eqs.(^) and (§) are 

n(t) = tz (9) 

x(t) = 0, 

representing a motionless particle at the origin with its magnetic moment 
(and spin) pointing antiparallel (n(t) = — z) to the direction of the field at 
that point and a similar solution but with the magnetic moment pointing 
parallel to the direction of the field (n(t) = +z). 
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3.2 Stability of the solutions. 

To check the stability of these solutions we now add first-order perturbations. 
We set 

n(t) = Tz+e x (t)x+e y (t)y (10) 
x{t) = + Sx(t), 

substitute these into Eqs.(^) and @, and retain only first-order terms. We 
find that the resulting equations for 5x(t), e x {t) and e y {t) are 

d 2 5x ^B' ± 

IF = ~ e " { > 

de y n 

We now look for oscillatory (stable) solutions for these equations. Setting 



inside Eqs.([TT|) yields 



= x e-^ 




= e x ,oe~^ 




— c p- iwt 

— £y,0 e 






( x 


»B /S ■ 




iuj J 


V e y,o 



u 2 

±HB'JS iu fxB /S j * j r,,,, J = { I . (12) 

—fiB Q /S iuo 




This equation has non-trivial solutions whenever the determinant of the ma- 
trix A vanishes. Thus, the secular equation 



*£detA = -W-^Y +(—\ '=Fl = (13) 

determines the eigenfrequencies u of the various possible modes. When the 
lower sign is taken in Eq . (|13"D , corresponding to a spin parallel to the mag- 
netic field, we find that one of the roots for uo 2 is purely negative. This 
indicates that one of the roots for u has a positive imaginary part for any 
K and hence, the solution is unstable. We concentrate now on the other 
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possible solution, corresponding to n(i) = — z: Here, the solutions for oj 2 are 
given by 



1 VI - 4:K 2 



2K 2 



± 



2K 2 



2 

u vib- 



(14) 



For K — > 0, the slow mode (minus sign) represents the vibration of the 
particle in the adiabatic potential, and the fast mode (plus sign) represents 
the precession of the spin about the z component of the magnetic field, as is 
shown explicitly by the form of the eigenvectors. The general form of these 
eigenvectors may be written in terms of an arbitrary amplitude parameter A 
as 



/ 5x \ 













B\ V lu J 

-1 



A 



(15) 



which for small K reduces to 




for the vibrational mode, and to 




prec 



( So \ 

-iK 
-1 



A 



(16) 



/ B T ^ 2 \ 



B\ 



-K 2 



V 



—I 

-1 



.4 



(17) 



/ 



for the precessional mode. From Eq.(|l^) we learn that in the vibrational 
mode, the amplitudes of the translational motion of the particle and the 
y-component of its spin are large compared to the amplitude of the x- 
component of the spin. Furthermore, the ratio (Sx /e y fi) vib = —B /B' ± shows 
that the amplitudes 5xo and e y $ are related in such a manner that the di- 
rection of the spin is antiparallel to the direction of the local magnetic field. 
Eq. (|i~7D for the precessional mode tells us that the amplitude of the transla- 
tional motion of the particle is negligible, thus in this mode the particle is 
essentially fixed and its spin precesses around the direction of the magnetic 
field at the origin. 
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Due to the coupling between the translational and the precessional de- 
grees of freedom, the mode frequencies given in Eq . (|T4]) , change with increas- 
ing K. 

A stable solution requires that all mode frequencies uj be real. Conse- 
quently, stability means that the roots for uj 2 are real and positive. From 
Eq.([L4]) we find that this happens when 

K < 0.5. (18) 

For the following discussion it is also useful to note that in the region K < 0.5, 

uj 2 ib < uj 2 < 2u vib ; for the vibrational mode (19) 
uj 2 > 2uj 2 ib ; for the precessional mode. 

Fig.0 shows the real and imaginary parts of the frequencies uj of the two 
modes as a function of K. We note that when K < 0.5, the imaginary parts 
of both frequencies vanish, indicating a stable solution. The mode with the 
lower frequency goes asymptotically to uj — > uj vib when K — > 0. This is the 
vibrational mode. The mode corresponding to the higher frequency, goes 
asymptotically to uj — ► uj prec = uj vib /K as K —* 0. This is the precessional 
mode whose frequency, as we mentioned earlier, is inversely proportional to 
the spin. When K > 0.5, both solutions for uj 2 possess an imaginary part 
and the system becomes unstable. 



3.3 The excitation energy of the modes. 

The excitation energy of a given mode is defined as the difference between 
the energy of the mode and the energy of the stationary state, 

1 / d5x x 2 



e = -^n-B+-m ^— j - (iB . (20) 

Note that the energy is bilinear in the coordinates and hence, one cannot 
neglect the z-component of the spin. Instead, one must set 



n ■ z 



Thus, the correct expression of the energy for small amplitudes is 



e - -V Q (4 + 4) B + B' ± Sxe y j ( 



d5x x 2 



(21) 
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positive imaginary part 
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Figure 1: Real and imaginary parts of the mode frequencies as a function of 
K. 
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In this expression, the modes have to be written in real form, 

Sx ( t ) = ^ (^.Y Acos{ut) , (22) 



e x (t) = -K—A sin (ut) , (23) 



(t) = -A cos (ut) . (24) 



Substituting Eqs.(H), (H) and (H) into Eq.® and using Eq.© one finds 
that 

6 _ 2 ^-^ 2 . (25) 



//-B 2cj 2 

Using Eq.(|T9D we conclude that for < K < 0.5, the excitation energy of the 
vibrational mode is positive while the excitation energy of the precessional 
mode is always negative. At the point K = 0.5, where the two modes coalesce, 
the excitation energy vanishes. We will further refer to these observations in 
the following section. 

3.4 The effect of viscous friction. 

When friction is introduced into the system, the equations of motion become 

d^x d dx 

m ^ = ^- B) - rt Tt (26) 



and 



rydh „ dh 

S— = un x B— r„nx — , 27 
dt ^ p dt 



where r t and r p are translational and precessional friction coefficients, respec- 
tively. The second term on the right-hand side of Eq.(|27D is the spin-damping 
contributed by the change in the direction of the spin from n to n+<in. Since, 
by definition, n is a unit vector, dn is perpendicular to n. Thus, Q± = \dh/dt\ 
is the angular velocity associated with the change of n. Since the direction 
of fl± must be perpendicular to both dh and n we form the cross product 
Q± = fix (dh/dt) which incorporates both the correct value and the right 
direction. Multiplying Q± by r p yields the spin-damping term. 
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To first order in r r and r t the secular equation in this case is given by 

K^ n 2i^ n K^ + ^ (J^) 2 r t K* - iu n (Jjf) \ t j + l + iu> n K T f = 

(28) 

where we defined 

_ UJ 

u) n = 

to make the expression simple. Let co> n o be the eigenfrequencies uj n of the 
frictionless problem, given by Eq.(|13j). When adding a small friction to the 
problem, the eigenfrequencies will change by a small amount 8u n . We find 
an approximate expression for 5ui n by expanding Eq. fl2"8|) around uj nt0 to first 
order in Su n and making use of Eq. fllTf ) . This gives 

iK ( y< (2< - 1) + n (B /B' ± f 
& " " M I, <^ 

Eq. (|29|) has an interesting consequence: From Eq.(|l9D we find that the 
numerator in Eq.(p9|) is positive for both modes while the denominator is neg- 
ative for the vibrational mode and positive for the precessional mode. We 
therefore conclude that friction, either translational or precessional, stabilizes 
the vibrational motion and, simultaneously, destabilizes the precessional mo- 
tion. The system all together becomes of course, unstable. 

The fact that spin damping leads to an exponential growth of the fast 
mode is no surprise in view of its negative excitation energy. Also, the expo- 
nential decay of the slow mode due to translational friction is to be expected 
on account of its positive excitation energy. What is important is the fact 
that due to the coupling between translation and precession, translational 
friction causes an exponential growth of the fast mode, with a growth time 
which, compared to the effect of spin damping, is smaller by a factor of 
r t K 2 S 2 / /j,mr p BQ in the limit of small K. 




4 Quantum-mechanical analysis. 

4.1 The Hamiltonian and its diagonalized form. 

In this section we consider the problem of a neutral particle with spin S = h/2 
in a ID inhomogeneous magnetic field from a quantum-mechanical point of 
view. Unlike the classical analysis, in which the derivation was valid for any 
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value of the adiabaticity parameter K, we concentrate here on the behavior 
of the system when K is small. We choose to analyze the case of a spin 
1/2 particle because this case already shows the essentials of the quantum- 
mechanical problem. 

Now, it is convenient to express the dependence of the magnetic field on 
x in terms of its amplitude B{x) and its direction 6{x) with respect to the z 
axis. Thus, Eq.@ is rewritten as 



B = B (x) (sin [9 (x)) y + cos [9 (x)) z) 



where 



B{x) 
9(x) 



B\ 1 



arctan 



B' ± x 
B' ± x 



(30) 



(31) 



The time-independent Schrodinger equation for this system is 

-h 2 d 2 



— fiB(x) (sin [9 (x)] a y + cos [9 (x)] a z 



2m dx 2 

where a y and a z are the Pauli matrices given by 



* (x) = E*f! (x) 



(32) 



0V 



-i 

1 



1 
-1 



E is the eigenenergy, and \P is the two-component spinor 

* = 

In matrix form Eq . (|32"D becomes 



■0T ( x ) 
*I>1 ( x ) 



(H K + H M ) 
where H K and H M , given by 

-h 2 



■0T fa) 
fa) 



(33) 



(34) 



\ 



2m 



H M = jj,B (x) 





dx 2 I 
-cos[0(a;)] isin[9(x)] 
-i sin [6 1 (x)] cos [6* (x)] 



(35) 
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are the kinetic part and the magnetic part of the Hamiltonian if, respectively. 

In order to diagonalize the magnetic part of the Hamiltonian, we make a 
local passive transformation of coordinates on the wave function such that the 
spinor is expressed in a new coordinate system whose z axis coincides with the 
direction of the magnetic field at the point x. We denote by R (x) the required 
transformation and set ty' = R^. Thus, 1 J r ' represent the same direction of 
the spin as before the transformation but using the new coordinate system. 
The Hamiltonian in this newly defined system is clearly given by RHR^ 1 . 
In the case of the magnetic field given in Eq . (|30|) , the required operation is 
a rotation by an angle 9 (x) around the — x axis. The operator that affects 
the wave function in this manner is |16| 



R 



exp 



e A 



while its inverse is given by 



R- 1 



exp 



' 9 ' 



cos (9/2) -i& x sm (9/2) 



cos (9/2) + icr x sin (9/2). 



It is easily verified that the transformation indeed diagonalizes the magnetic 
part of the Hamiltonian, 



RHmR 1 
For the kinetic part we find 

d 2 1 



-/i-B (x) a z . 



RHkR 1 



2m 



dx 2 



dOy JdBd 1 d 2 9 
A\dx) \dxdx 2dx 2 



Or 



, .(dOd ld 2 9\ A TT . . 

Note that i — I r^7 a x IS Hermitian. 

\dx ox 2 dx 2 J 

Thus, the Hamiltonian of the system in the rotated frame may be written 



as 



where 



H = H, 



diag 



int 



(36) 



Hdiag 



2m 



d 2 _ 1 fd9 

dx 2 4 \dx 



hit 



fiB (x) a z 



(37) 



. h 2 fd9 d ld 2 9" 
2m \dxdx 2 dx 2 ; 



13 



The first part of the Hamiltonian Hdiag is diagonal. It contains the kinetic 
part ~ d 2 /dx 2 , a term whose form is =F fiB(x) which is to be identified as 
the adiabatic effective potential, and a term ~ (d6/dx) 2 which appear due 
to the rotation. The second part of the Hamiltonian H int contains only non- 
diagonal components. These will be shown to be of order O (K) and hence 
may be regarded as a small perturbation. We proceed to find the eigenstates 

Of Hdiag- 



4.2 Stationary states of Hdiag- 

Since Hdiag is diagonal, the two spin states of the wavefunction are decoupled. 
We then seek a solution of the form 



^=UwJ '• E=El ' 

referred to as the spin-down state, and another solution 

*t=( ^i x) \ ; E — Et 







(3* 



(39) 



which we call the spin-up state. 

The equation for the non-vanishing component of the spin-down state is 
given by 



h 2 ( d 2 l fdd\ 2 . 



V>1 = E i^h 



(40) 



whereas the equation for the non-vanishing component of the spin-up state 
is 



2 



de\ 

2m \ dx 2 A\dx J 



— fiB (x) 



= E-^ipp 



(41) 



We now show that in the limit of small K we can neglect the term ~ 
(dO/dxy in both Eq.(gOD and Eq.(|l]): We compare the order of magnitude 
of the term fiB (x) to that of the term h 2 (d8/dx) 2 /8m. Using Eq . (|3l"l) it 
can be easily shown that the maximum value of d9/dx is B'j_/ Bq whereas the 
minimum value of fiB (x) is (iBq. Thus, 



fiB (x) 



h 2 fdex 

8m \dx J ] 



8fimBl 2 
(B' L ) 2 K 2 ~ K 2 ' 
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and we reach the conclusion that when K is small enough we can neglect the 
term ~ (d6/dx) 2 . Under this approximation, Eqs. (|40|) and ( |41~D simplify to 



and 



2m dx 2 



tf_cP_ 
2m dx 2 



+ \xB [x] 



- fiB (x) 



■01 = E^ipp 



(42) 



(43) 



The approximate solutions of these equations are outlined in the next two 
subsections. 



4.2.1 Stationary spin-down states. 

Eq.(^) represents a particle in a symmetric attractive potential. If the extent 
of the wave function is small enough, we can expand B (x) to second order 
in x 



B (x) ~ B 



2 I B 



(44) 



and apply the well-known solution of a harmonic oscillator in one dimension. 
We now derive the condition for which this approximation is valid: We recall 
that the ground-state wave function of the harmonic oscillator is given by||17|| 



2mio vi } ) x 2 



(45) 



The extent of this wave function over which it changes appreciably is given 
by 



Ax i 



h 



2mu vib 

whereas the extent over which fiB (x) changes significantly (see Eq. ( 

Ax mB ~ B /B' ± . 
Thus, the ratio between these two length scales is 

Ax i 



Ax 



fj,B 



2mu 2 ib B 2 



K. 



(46) 

-])) ^ 

(47) 

(48) 
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We therefore conclude that when K is small enough, the harmonic approxi- 
mation is justified. The wave function ip^ (x), given by Eq . (fl5"l) , then repre- 
sents the lowest possible bound state for this system. This state corresponds 
to a trapped particle. The energy of this state is clearly 

h 

E i = fiB + -u mb = fiB (1+K)~ fiB , (49) 
while its full spinor representation is 



/ \ 

2muj V it ) x 2 
(mLo vib \V* r 



(50) 



4.2.2 Stationary spin-up states. 

Eq.(|^) describes a particle in a repulsive potential. It corresponds to an 
unbounded state representing an untrapped particle. In this case there is a 
continuum of states, each with its own energy. As we are interested in non- 
radiative decay, we focus on finding a solution with an energy which is equal 
to the energy found for the trapped state, that is 

£ T = E t ~ fiB . (51) 

Since H diag is symmetric under x — > —x, to each energy in the continuous 
spectrum there belongs one state with even parity and one state with odd 
parity. Further, H int connects states with opposite parity only. Therefore, 
since the spin-down ground state is even in x, we need the spin-down state 
with odd parity. 

When evaluating the lifetime in the next section, we compute the matrix 
element of H int between the states ip-\(x) and ip[(x). Thus, most of the contri- 
bution to this integral comes from the region in x where ipi{x) is substantial. 
According to Eq.(f48|), fiB(x) changes very little in this range and, as a first 
approximation, we may take the potential in this region as uniform, 

/iB(x) ~ fiB (52) 

in Eq.fl43D, and then its solution becomes 

r h { x) = Cmn(^^xY (53) 
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with the full spinor representation 

= ( Can {~^~ X J | (54) 

The approximate wave function given in Eq.([5tJ) is periodic near the ori- 
gin, and has a period of 

Ax T = -===. (55) 
Comparing it to Ax± given in Eq. (|46|) , we find that 



Ax i 



if, (56) 



which shows that, when if C 1, the wavefunction if)^ (x) executes many 
oscillations in the region where ip^ (x) is appreciable. 

4.3 The lifetime. 

To evaluate the lifetime T esc of the particle in its trapped state, which is the 
average time it takes for the particle to escape, we calculate the transition 
rate from the bound state given by Eq. (|50D , to the unbounded state Eq. fl54|) , 
according to Fermi's golden rule|L9]]. Thus, 



1 2tt , , , |2 



Tesc 

where 



tH^wI^) (57) 



+oo +oo 

H U = = jd^l (,) g) (ig + ig|) * (,) 

— oo — oo 

(58) 

is the matrix element of Hi nt Eq.(^) between and and g(E^) is the 
density of the final states at energy E^. 
Using Eq . (|53|) and Eq . (|31~D we find that 

id^ j (SlY 



dOdijj^ B' ± (y/^m\iB K 



dx dx B \ h 
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Thus, when K is small we may neglect the contribution of d 2 6/dx 2 to the 
integral Eq . fl58|) and then we find that 



H 



2m 



+00 



dxipl (x) 



d9 dt/ji (x) 
dx dx 



(59) 



The integrand in Eq. fl59|) consists of a product of three functions: The 
function ip* whose 'width' is about Ax± (given in Eq.(f46])) around the origin, 

the function d6/dx whose extent around the origin Ax^b is roughly ^Jl/K 
larger than Ax± and the function difi^/dx which is a periodic function with 
a characteristic period Axj which is \JK smaller than Ax±. This suggests 
that we can approximate the integral in Eq . (|59|) by substituting d9/dx for 
its value at x = 0, 



d6 
dx 



El 

B 



Substituting Eqs.©, © and (||) into Eq.© gives 



H 



h c /2/i(Si 



where we have used the definite integral 



1/4 



exp 



2hui. 



vib 



+00 



cos (bx) dx 



exp 



4a 



(60) 



(61) 



(62) 



When Eq.© is substituted into Eq.fl57|) the term \C\ g(E^) appears. 
This term can be calculated by temporarily introducing suitable boundary 
conditions: Assume that the system is bounded by an infinite potential wall 
at x = ±L/2, the length L being large compared to Ax^ yet small when 
compared to Ax^b- In this case, the uniform potential approximation still 
holds, and we may use the well-known result for the density of states for a 
particle in a ID infinite potential well, 



g(Ef = fiB ) 



mL 2 



AfiB 7T 2 h 2 ' 

The evaluation of the normalization constant C gives 

\C\ 2 = 2/L, 



(63) 



(64) 



18 



and therefore 



\C\ 2 g(E^ 



m 



flB 7T 2 h 2 ' 

Finally, using Eqs.(^) and fl6"l|) inside Eq. (|57|) gives 



(65) 



u vib V 2nK 



exp 



2K 



T, 



nib 



(2ti) 3 K 



exp 



2K 



(66) 



where T vi b = is the period of classical oscillations inside the trap. 

Looking back at the calculation of the matrix element, Eq. ([[]]), one re- 
alizes that for small K the wave function ip^ oscillates strongly in the region 
where the wave function tpi is essentially different from zero. Thus, successive 
subintervals compensate each other very effectively, and it is possible that 
the regions \x\ > Ax^ where the x-dependence of the potential becomes im- 
portant, contribute more effectively to the integral than the simple estimates 
indicate. We have therefore carried out a more accurate calculation, using 



the WKB approximation^^], which takes the x-dependence of the potential 
into account. The calculation, which is given in the Appendix, yields the 
result 



J- esc — 



1 33 



uj v ib 32 V 2ttK 



exp 



2 fl 
— arctan - 
K V 4 



(67) 



which does indeed differ from Eq.(^) not only by the pref actor, being unity 
in Eq.(|66"D and 33/32 in Eq. (|6T|) , but even by the exponent which is 0.5 in 
Eq.® and 2 arctan (1/4) ~ 0.489 96 in Eq.(|7|). 



5 Discussion. 

Summarizing all we have found we conclude that the problem we have studied 
has three important time scales: The shortest time scale is T prec , which is 
the time required for one precession of the spin around the axis of the local 
magnetic field. The intermediate time scale is T vi b = T prec /K, which is the 
time required to complete one cycle of the center of mass around the center 
of the trap. These two time scales appear both in the classical and the 
quantum-mechanical analysis. The longest time scale (provided K is small) 
T esc , which is not present in the classical problem, is the time it takes for the 
particle to escape from the trap. 

Whereas the classical analysis yields an upper bound of K = 0.5 for 
trapping to occur, no such sharp bound exists in the quantum-mechanical 



19 



analysis. This is related to the fact that one cannot associate an effective po- 
tential well with a finite barrier with the system. Nevertheless it is interesting 
to compare the classical bound with the value of K for which the exponent 
in the expression for the quantum- mechanical lifetime becomes equal to 1: 
For the uniform- field approximation, Eq . (|66|) , the two values happen to agree 
exactly, and in the WKB approximation, Eq. (|67D , they differ only by a small 
amount. Thus, the quantum- mechanical condition for trapping to occur is 
essentially the same as the classical condition. 

As an example, we apply our results to the case of a neutron and an atom 
trapped with a field B Q = 100 Oe and B Q /B' X = 10cm. These parameters cor- 
respond to typical traps used in Bose-Einstein condensation experiments [|7|]. 
The results, being correct to within an order of magnitude, are outlined in 
the following table: 





B = 100 Oe 




B /B' ± 


= 10cm 




Neutron 


Atom 


m gr 


~ io- 25 


~ IO" 22 


H emu 


~ io- 23 


~ IO" 20 


K 


~ 10~ 5 


~ io- 8 


Tpyec S6C 


~ IO" 6 


~ io- 9 


Tyib SGC 


~ 10" 1 


~ 10" 1 




~ io( 104 ) 


~ io( 107 ) 



We note that in both cases K is very small compared to 1. Consequently, the 
calculated lifetime of the particle in the trap is extremely large, suggesting 
that the particle (either neutron or atom) is tightly trapped in this field. 

Though the toy model presented in this paper is very simple, preliminary 
studies |y| show that behavior similar to what we found in this model trap 



appears also in more realistic magnetic traps such as those used in Bose- 
Einstein condensation experiments. For example, a similar analysis of a two- 
dimensional Ioffee trap|2j shows similar behavior when the classical analysis 
and quantum-mechanical analysis are compared, but in this case we find that 
T esc ~exp[2/^ 3 ]. 

The problem studied in this paper deals with a spin 1/2 particle. Though 
this fact has little influence on the solution of the classical problem, the 
extension to higher spin values complicates the analysis of the quantum- 
mechanical problem. In this case one has to deal with a (2S + l)-component 
spinor, and the interaction Hamiltonian does no longer connect the (-S)- 
state to the (+S')-state, but only to the (-S + 1) and (— S + 2) states which 
for S > 5/2 will still be trapped. 
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A Calculation of ip^ and T esc by the WKB 
method. 



In this appendix, we outline the solution of Eq.(|43|) by the WKB approx- 
imation and calculate the resulting lifetime. The validity of the WKB ap- 
proximation is guaranteed by the following argument: We first note that, for 
Ei = /i-Bo; the local de-Broglie wavelength A (x) is given by 



A(x) 



2iih 



P E A X >\&E-r=flB 



(6* 



where 



Pe } {x) = J2m(E^ + /jB (x)) 



2mfiB 2 + 



1 (B\x 



2 V B 



(69) 



is the classical momentum of the particle. Thus, the rate of change of A over 
x is 



dX 



dx 



2irhfim [2m(fiBo + fiB (x) 



-3/2 



dB(x) 



dx 



(70) 



Its maximum value is given by 



dX 



dx 



\ivah 



[2m({iB + fiB )} 



3/2 



/jtuBq 



(71) 



Hence, in the adiabatic approximation, where K <§C 1, the validity condition 
of the WKB approximation \dX/dx\ 1 is satisfied. 
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In this approximation the wavefunction ipi(x) corresponding to the energy 
is given by []T8 



D 



Pe 1 W 



sin | p E] .(x')dx' J , 



(72) 



where a, the lower limit of the integration, is yet undetermined. 

We still need to evaluate the product \D\ 2 g(E^ = fiB) which appears 
when calculating T esc by Eq.(^). We therefore temporarily introduce bound- 
ary conditions by assuming that the system is bounded by an infinite poten- 
tial wall at x = ±L/2. Thus, by setting a = —L/2 the boundary condition 
at x = —L/2 is automatically satisfied. We note that we are looking for a 
highly excited state ipi an d we consider the eigenstate corresponding to the 
nearest higher energy E^ + AE. The latter must fulfill the requirement that 
the phase of ^ at x = L/2 changes by it when going from E^ to E<\ + AE. 
Thus, 



1 

h 



L/2 



L/2 



PE r +AE(x')dx' 



L/2 



L/2 



p E . (x')dx' 



7T. 



(73) 



Since E-\ 3> iiB{x) everywhere we may write the difference in Eq.(|T3"D as a 
derivative 



h dE r 



L/2 



L/2 



p E ^(x')dx' = 7T. 



Using Eq.(^) we rewrite Eq.([73) as 



\ae 
h 



L/2 



m 



-dx' 



7T. 



-L/2 PE^X') 

Normalization of the wavefunction Eq. (|72|) , gives on the other-hand 



L/2 



L/2 



|-?/>f (x)| 2 dx ~ \D\ 



L/2 



1 



L/2 \PE r (x) 
L/2 1 



sm 



/ 

1 

h 

\ -L/2 



\ 



p E (x')dx' 



J 



dx = 1, 



-L/2 PE{X) 



(74) 



(75) 



dx (76) 



where we have neglected the terms containing fast oscillations by replacing 
the sin 2 () term by its average value 1/2 . Substituting Eq.flTB]) into Eq.(|75|) 
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gives 



, ,9 1 , ,9 , 2m 
\D\ — — = \D\ g(E-t) = — 
1 1 AE 1 1 yv u iih 



(77) 



where S'(-E'f) = 1/AE is the density of states at energy Ep 

To evaluate the integral in Eq. (|59|) we use the stationary phase method|20 
Using Eqs.([72|), (|45|) and ( |60"D we rewrite this integral as 



where we have defined 



-too 

I± = A J dxg(x) exp/ ± (x) 



~ 4m\ 27vh J V ^ / V 



and 



= 



/±C«0 



2ma; t; jbX 2 
ft 



i 7 I (a/) da/. 



(75 



Note that the lower integration limit in /± (x) has been chosen as to make 
the function ^ (x) antisymmetric under x — > —x. 

We continue to work out I + first: According to the stationary phase 
approximation we should first find the point x s for which the phase f+(x) is 
stationary. This is accomplished by solving 



d 



dx dx 



ImiOvibX 1 

h 



+ ^ / PE r (a/) dx' 



(79) 



for x. Using Eqs.(51) and ( fHj) and after a few algebraic steps we find that 
the only solution to Eq.(|79f) is 



x: 



17 B\ 



*0) 
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Now we approximate the integral 1+ as 



+00 



A \ dxg(x) exp [f+(x)\ ~ Ag(xj), 



2tt 



exp [A 



The evaluation of f + (xf) is performed in the following way: 



IX' 



iW4/17B /B' ± 



2mu vib (x+ 

h 



1 

h 



p E ^ (x) dx 



2 y/4/17B /B' ± 
2 



17 B\ 



h 



Pe^ (iy) dy 



ft V m5 i V 17 5', 







h 



2mfiB 2 - 



1 (B\y 



2 V fin 



arcsm 



1 



(B' ± yh 2 \vn 

where we used the definite integral 



arctan (1/4) 
K 



/ \J a 2 — x 2 dx = ~y\j a 2 — y 2 + -a 2 arcsin (—~\ ; for y < a. 
Jo 2 2 \a/ 

The term f+(xf) is calculated straightforward with the result that 

1 



m\i (B' ± 
B 



-2-i 



17y/l ~ 1/17 



(82) 



m/i (B 



I \2 



1089 



exp 



1 arctan 



17 



136 



B V 1088 

Substituting this together with g(xf) and f(xf) into Eq.([H]) gives 
B (17) 1/4 [to 



Ba 



B' L 2 V \7\\mix{B 



I) 2 



(83) 



x 



1088 
1089 



1/4 



exp 



17 



-1 arctan Q) -^arctan ^ 
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Likewise, a similar calculation for J_ gives 



/_(*7) = /+(*+) 
#00 = -fl'(^); 



leading to the result that 

i_ = -%A 
1088 



1/4 



o 



X 



Bp (17) 

1/4 




1089 



cxp 



— arctan - 
K V4 



arctan 



By repeating the steps that led to Eq.flBq) we find that 

1 33 



1 



32 V 2nK 



cxp 



" 2 / 




— arctan 


'M 


K 
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